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Abstract — In  this  paper,  we  address  the  deployment  problem 
for  differentiated  detection  requirements,  in  which  the  required 
detection  probability  thresholds  at  different  locations  are  dif¬ 
ferent.  We  focus  on  differentiated  deployment  algorithms  that 
are  applied  to  the  probabilistic  detection  model,  since  it  is 
more  realistic  than  the  binary  detection  model.  We  show  that 
the  relationship  between  the  node  deployment  strategy  and  the 
logarithmic  collective  miss  probability  distribution  is  Linear  Shift 
Invariant  (LSI).  Using  this  property,  we  formulate  the  differen¬ 
tiated  deployment  problem  as  an  integer  linear  programming 
problem,  which  is  a  well  known  NP-hard  problem.  We  propose  a 
differentiated  node  deployment  algorithm  called  DIFF_DEPLOY, 
which  achieves  much  better  performance  than  the  state-of-the-art 
node  deployment  algorithm  for  both  uniform  and  differentiated 
detection  requirements. 

I.  Introduction 

The  emergence  of  wireless  sensor  networks  gives  rise  to 
many  applications,  such  as  military  surveillance  [1]  and  habi¬ 
tat  monitoring  [2].  In  these  applications,  it  is  an  important 
requirement  to  provide  adequate  sensing  coverage  to  achieve 
acceptable  quality  of  service. 

In  previous  papers  on  sensing  coverage  [3][4][5][6],  a  binary 
detection  model  is  assumed.  In  a  binary  detection  model, 
sensor  node  can  detect  a  target  with  a  100%  probability 
provided  that  the  target  is  within  its  sensing  range,  and  cannot 
detect  a  target  beyond  the  range.  In  this  paper,  we  use  a 
probabilistic  detection  model,  because  in  reality  the  detection 
of  a  target  is  not  deterministic  due  to  the  uncertainty  associated 
with  sensor  detections.  When  a  target  is  within  a  sensor’s 
sensing  range,  whether  the  target  is  detected  by  the  sensor 
is  probabilistic. 

With  a  probabilistic  detection  model,  we  can  hardly  have 
a  100%  detection  probability  for  all  the  geographic  points  in 
the  target  area.  In  many  surveillance  applications,  the  system 
requires  different  degrees  of  surveillance  at  different  locations. 
The  system  may  require  extremely  high  detection  probabilities 
at  certain  sensitive  areas.  However,  for  some  not  so  sensitive 
areas,  relatively  low  detection  probabilities  are  required  to 
reduce  the  number  of  sensors  deployed  so  as  to  decrease  the 
cost.  In  this  case,  different  areas  require  different  densities  of 
deployed  nodes.  When  the  detection  probability  thresholds  at 
different  subareas  are  specified,  the  minimum  number  and  the 
deployment  locations  of  the  sensors  need  to  be  decided  based 
on  the  specification  and  the  probabilistic  detection  model.  We 
name  it  as  the  “Differentiated  Deployment”  problem.  The  goal 


of  the  paper  is  to  develop  a  deployment  strategy  to  satisfy  the 
different  detection  probability  thresholds  at  different  locations 
using  minimum  number  of  nodes. 

Minimizing  the  number  of  nodes  deployed  may  not  be 
critical  when  the  cost  of  node  is  negligible  and  we  have 
excessive  number  of  nodes.  However,  as  long  as  the  prices  of 
sensor  nodes  are  not  negligible,  to  reduce  the  number  of  nodes 
deployed  is  always  necessary.  For  example,  while  MicaZ 
motes,  which  are  commonly  used  in  terrestrial  environments, 
cost  over  120  dollars  each  [7],  some  sensors  used  in  undersea 
surveillance  applications  are  tens  of  or  even  hundreds  of  times 
more  expensive.  While  it  is  expected  that  the  cost  of  sensors 
will  decrease  as  the  technologies  advance,  it  may  not  happen 
quickly,  since  the  advances  of  the  technologies  also  result  in 
more  powerful  features  being  incorporated  into  a  single  node. 
This  is  especially  true  for  undersea  sensor  nodes.  Another 
issue  is  stealthiness.  By  minimizing  the  number  of  sensors 
deployed,  the  risk  of  the  system  being  detected  by  adversaries 
is  also  reduced.  Even  if  the  cost  and  the  stealthiness  are  not  an 
issue,  knowing  the  minimal  number  and  the  deployment  loca¬ 
tions  of  the  nodes  provides  us  the  guidance  on  how  to  deploy 
the  redundant  sensors  to  improve  reliability.  For  example,  the 
redundant  sensors  can  be  deployed  proportionally,  based  on 
the  density  of  the  sensors  after  a  minimum  number  of  sensors 
are  deployed. 

We  assume  that  we  have  good  control  of  the  node  deploy¬ 
ment,  i.e.,  we  can  place  the  sensors  into  the  exact  targeted 
locations,  either  manually  [8]  or  air-dropped  [9].  However, 
when  spatial  error  is  inevitable  in  the  node  deployment,  we 
can  model  the  uncertainty  of  nodes’  actual  locations  using 
Gaussian  distribution  [10].  We  show  that  to  consider  uncer¬ 
tainty  in  the  node  deployment,  we  only  need  to  modify  the 
node  detection  model.  Therefore,  our  algorithm  is  applicable 
to  both  precise  node  deployment  and  node  deployment  with 
uncertainty. 

In  this  paper,  we  show  that  the  relationship  between  the 
node  deployment  strategy  and  the  logarithmic  collective  miss 
probability  distribution  of  the  sensor  field  is  Linear  Shift 
Invariant  (LSI).  By  taking  advantage  of  this  property,  we 
formulate  the  node  deployment  problem  as  an  integer  lin¬ 
ear  programming  problem,  which  is  a  well  known  NP-hard 
problem.  Further,  based  on  the  linear  relationship,  we  devise 
a  differentiated  deployment  algorithm  called  DIFF  DEPLOY, 
which  outperforms  the  state-of-the-art  probability  based  node 


1-4244-0626-9/06/$  20.00  (c)  2006  IEEE 

This  full  text  paper  was  peer  reviewed  at  the  direction  of  IEEE  Communications  Society  subject  matter  experts  for  publication  in  the  IEEE  SECON  2006  proceedings. 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

2QQ£  2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2006  to  00-00-2006 

4.  TITLE  AND  SUBTITLE 

Deployment  Strategies  for  Differentiated  Detection  in  Wirelesss  Sensor 
Networks 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

University  of  Virginia, Department  of  Computer  Science, 151  Engineer’s 
Way, Charlottesville, VA, 22904-4740 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

18.  NUMBER  19a.  NAME  OF 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE 

unclassified  unclassified  unclassified 

10 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


deployment  algorithm  MIN_MISS  for  both  uniform  and  dif¬ 
ferentiated  detection  requirements. 

The  remainder  of  the  paper  is  organized  as  follows.  We 
review  related  work  in  Section  II.  We  describe  the  sensor 
and  terrain  models  in  Section  III.  Section  IV  presents  how 
to  relate  the  detection  probability  to  system  performance.  We 
formulate  the  node  deployment  problem  in  Section  V.  In  Sec¬ 
tion  VI,  we  present  our  differentiated  deployment  algorithm 
DIFF_DEPLOY.  Then  we  briefly  describe  MIN_MISS  and 
the  way  uncertainty  is  handled  in  Section  VII.  Section  VIII 
contains  a  complete  evaluation  of  different  deployment  al¬ 
gorithms.  We  present  the  conclusions  and  future  work  in 
Section  IX. 

II.  Related  Work 

In  sensor  network  based  surveillance  applications,  it  is 
critical  to  determine  where  to  deploy  the  sensors  to  satisfy 
the  sensing  coverage  requirement.  A  large  number  of  pub¬ 
lications  exist  in  the  literature  [  10] [3] [  1 1  ] [4] [5] [  12] [6] [13], 
However,  most  of  them  are  based  on  binary  detection  models 
[3] [4]  [5] [6], 

When  the  binary  detection  model  is  used,  if  we  only 
consider  sensing  coverage,  to  determine  the  minimum  number 
of  nodes  to  cover  the  whole  area  is  actually  a  circle  covering 
problem  [5][14],  A  lot  of  work  also  considers  communication 
connectivity  requirements  while  minimizing  the  number  of 
nodes  deployed  [5] [6]. 

Although  the  binary  detection  model  simplifies  the  analysis, 
it  is  not  realistic  in  many  cases.  In  reality,  a  sensor  detection  is 
a  probabilistic  event,  in  which  a  target  is  not  always  detected 
by  a  sensor  due  to  its  limited  sensing  capability  and  processing 
power.  In  [12],  Meguerdichian  et.  al.  assume  that  the  sensor 
coverage  decreases  as  the  distance  from  the  sensor  increases 
and  propose  an  algorithm  based  on  Voronoi  diagram  and 
Delaunay  triangulation  to  find  the  maximal  breach  path  and 
maximal  support  path.  Although  that  algorithm  can  be  used  as 
a  guide  to  deploy  extra  sensors  to  improve  breach  coverage, 
it  can  not  be  used  for  differentiated  deployment  problem 
unless  both  different  node  detection  models  and  differentiated 
detection  requirements  are  considered.  In  [11],  Dhillon  et. 
al.  propose  to  use  probabilistic  detection  models  to  better 
address  this  problem.  In  these  models,  the  probability  of 
target  detection  by  a  sensor  decreases  exponentially  when  the 
distance  between  them  increases. 

Several  deployment  algorithms  [10][1 1][  13]  are  proposed 
when  the  probabilistic  detection  models  are  used.  The  algo¬ 
rithm  described  in  [13]  assumes  the  mobility  of  the  sensors, 
while  we  consider  static  sensors  that  do  not  move  once  they  are 
deployed.  The  deployment  algorithms  described  in  [  10]  [  1 1  ] 
are  similar,  between  which  MIN'  MISS  proposed  in  [10]  is 
more  sophisticated.  MIN_MISS  is  mainly  used  for  uniform 
detection  requirements,  in  which  the  specified  detection  prob¬ 
ability  threshold  for  the  entire  surveillance  area  is  the  same. 
Different  from  MIN_MISS,  DIFF .DEPLOY  is  designed  for 
both  uniform  and  differentiated  detection  requirements.  Fur¬ 
thermore,  Dil  l-  DEPLOY  achieves  much  better  performance 


in  both  cases. 

III.  Sensor  and  Terrain  Models 

In  this  paper,  we  use  a  two-dimensional  grid  to  cover  the 
sensor  field.  We  only  consider  the  detection  probability  at  the 
grid  points.  The  granularity  of  the  grid  is  adjusted  based  on 
the  precision  we  want  and  the  time  and  space  we  can  afford 
for  computation.  Suppose  the  dimension  of  the  grid  is  U  x  V. 
The  deployment  strategy  of  the  sensor  field  can  be  represented 
by  a  U  x  V  matrix  D,  in  which  D(x,  y )  denotes  the  number  of 
nodes  deployed  at  grid  point  (x,y).  If  D(x,y)  =  1,  it  means 
that  one  sensor  is  deployed  at  grid  point  ( x,y ).  If  D(x,y)  = 
0,  it  means  no  sensor  is  deployed  at  that  location.  Typically, 
D(x,y )  is  either  0  or  1.  However,  in  very  rare  cases  when 
the  detection  requirement  is  so  high  that  more  than  one  node 
need  to  be  deployed  in  a  small  neighborhood,  it  is  possible  that 
D(x,y)  takes  a  value  greater  than  1.  We  use  d((?n,?i), 
to  denote  the  Euclidean  distance  between  locations  (?n,  n)  and 

(m)- 

In  reality,  it  is  challenging  to  obtain  an  exact  and  precise 
node  detection  model,  because  the  detection  of  a  target  de¬ 
pends  on  many  factors.  We  classify  the  factors  into  three 
categories.  Factors  in  the  first  category  are  those  related  to 
the  sensors,  such  as  the  type  and  the  quality  of  the  sensor, 
and  the  sensing  algorithm  used  to  detect  the  target.  Factors 
in  the  second  category  are  those  related  to  the  target.  They 
include:  the  type,  size  and  shape  of  the  target;  the  duration 
of  a  target  staying  at  a  certain  location,  which  is  related  to 
the  speed  of  the  target;  the  distance  between  the  target  and 
the  sensor.  Factors  in  the  third  category  are  those  related  to 
the  environment  where  sensors  are  deployed.  This  category 
includes  the  following  factors:  the  surrounding  obstacles;  the 
weather,  such  as  the  temperature,  humidity  and  whether  it  is 
rainy  or  windy;  background  noise,  such  as  the  magnetic  field  of 
the  earth.  Almost  all  these  factors  are  experienced  in  real  field 
experiments  [15].  Many  of  these  factors  are  hard  to  model. 
Therefore,  it  is  very  difficult  to  get  a  totally  accurate  detection 
model. 

However,  it  is  possible  to  obtain  a  conservative  detection 
model,  as  what  has  been  done  to  obtain  a  conservative  sensing 
range  [16].  This  kind  of  model  can  be  obtained  by  using 
conservative  values  in  theoretical  analysis  or  through  extensive 
experiments.  For  example,  we  can  run  the  experiments  many 
times  under  different  weather  conditions,  and  use  the  lowest 
detection  probability  if  the  results  under  different  weather 
conditions  are  different.  In  this  way,  we  rule  out  several 
hard  to  model  factors  and  the  obtained  model  provides  the 
worst  case  detection  probability.  Although  it  may  be  time 
consuming,  it  is  possible  to  get  a  conservative  detection  model 
with  current  technologies.  Obtaining  the  detection  model  is 
not  only  useful  for  the  node  deployment  problem,  but  helps 
us  better  understand  the  performance  of  the  system.  Providing 
a  totally  accurate  detection  model  is  out  of  the  scope  of  this 
paper. 

In  [10]  [13],  a  probabilistic  node  detection  model  consider¬ 
ing  the  factor  of  the  distance  between  the  target  and  the  sensor 
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is  used,  which  is  shown  in  Equation  (1): 

p((m,  n ),  (i,j))  = 

e-arf((m,n),(i,j))  if  d((m,  n),  (i,  j))  <=  rs  ... 

0  if  d((m,n),(i,j))>rs 

In  the  above  equation,  p((m,  n),  ( i,j ))  is  the  probability  of 
a  target  at  grid  point  (m,  n)  being  detected  by  the  node  at  grid 
point  ( i,j),rs  is  the  sensing  range,  and  a  is  a  parameter  related 
to  the  physical  characteristics  of  the  sensing  device,  which  can 
be  obtained  from  field  experiments.  This  model  assumes  that 
the  detection  probability  varies  exponentially  with  the  distance 
between  the  target  and  the  sensor.  Although  this  model  is 
not  accurate,  it  reasonably  characterizes  the  behavior  of  range 
sensing  devices,  such  as  infrared  and  ultrasound  sensors  [17], 
This  equation  does  not  take  into  consideration  of  the  time 
duration  a  target  stays  at  a  certain  location.  In  reality,  if  a 
target  stays  at  a  certain  location  for  a  longer  time,  it  results  in  a 
higher  probability  of  being  detected.  Therefore,  we  incorporate 
the  factor  of  the  duration  into  the  node  detection  model: 

p((m,n),  ( i,j),t )  = 

1  -  I1  -  P((m.  n),  l))rt/T1  ift>T  , 

p((m,  n),  ( i ,  j))  if  0  <  t  ^  T 

In  this  model,  we  use  discrete  time.  The  time  unit  T  is  the 
time  period  during  which  the  sensing  algorithm  is  executed 
and  a  decision  is  made  on  whether  a  target  is  detected,  t 
denotes  the  duration  the  target  stays  at  location  This 

model  assumes  that  the  target  detection  in  different  time  units 
is  independent.  We  also  assume  that  1)  = 

Note  that  the  detetion  probability  at  a  certain  location 
for  one  time  unit  can  be  easily  converted  to  the  detection 
probability  for  multiple  time  units.  If  the  detection  probability 
for  one  time  unit  is  w,  the  detection  probability  for  l  time 
units  is  1  —  (1  —  w)1.  For  simplicity,  we  assume  that  the 
detection  probability  thresholds  specified  by  the  users  are  the 
probabilities  for  one  time  unit.  To  be  consistent  with  the  detec¬ 
tion  probability  thresholds,  all  the  probabilities  mentioned  in 
the  following  sections  are  the  probabilities  for  one  time  unit. 
Since  p((m,n),(i,j))  =  p((m,n),(i,j),  1),  we  directly  use 
p((m,ri),  ( i,j ))  instead  of  p((m,ri),  ( i,j ),  1)  for  brevity. 

Similar  to  p((m,n),  (■ i,j )),  pmiss((m,n),  ( i,j ))  is  used  to 
denote  the  probability  of  a  target  at  grid  point  (m,  n)  being 
missed  by  the  node  at  grid  point  (i,j),  which  is  shown  in 
Equation  (3): 

Pmiss((m,  n),  ( i,j ))  =  1  -  p((m,n),  ( i,j ))  (3) 

We  use  a  U  x  V  matrix  M  to  denote  the  collective  miss 
probability  distribution  of  the  whole  field.  The  collective  miss 
probability  M(x,y)  means  the  probability  of  a  target  at  grid 
point  (x,  y)  being  missed  by  all  the  sensors  deployed  in  the 
field,  which  is  given  by: 


(i,j)€.Grid 

=  f[  (1-p((a;.y).(i.i)))D<,’5)  (4) 

Grid  in  Equation  (4)  denotes  the  whole  sensor  field  where 
we  want  to  deploy  the  sensors. 

IV.  Detection  Probability  vs  System  Performance 

In  the  differentiated  deployment  problem,  realizing  the  rela¬ 
tionship  between  the  detection  probability  and  the  system  per¬ 
formance  helps  the  designer  to  specify  the  detection  probabil¬ 
ity  thresholds  based  on  the  different  performance  requirements 
at  different  locations.  In  this  section,  we  briefly  discuss  how 
to  infer  the  system  performance  from  the  detection  probability 
and  vice  versa. 

In  surveillance  applications,  primary  performance  require¬ 
ments  include  detection  performance  and  tracking  perfor¬ 
mance.  The  detection  probability  is  closely  related  to  both  of 
them.  The  detection  performance  is  about  whether  a  target  can 
be  detected  when  it  appears  in  the  network,  and  the  detection 
delay  if  the  target  is  detected.  The  tracking  performance 
represents  how  well  the  target’s  trajectory  is  obtained,  given 
that  the  target  is  detected. 

For  different  target  types,  we  use  different  requirements. 
We  divide  the  targets  into  two  categories:  static  targets  and 
moving  targets.  For  the  first  category,  we  are  mainly  interested 
in  the  detection  performance.  The  detection  probability  at  a 
certain  location  indicates  whether  a  target  can  be  detected  if 
the  target  stays  there  for  a  certain  period  of  time.  We  use 
Pdetect{x,y)  to  denote  the  detection  probability  at  grid  point 
(x,  y)  when  the  target  stays  at  (. x ,  y)  for  one  time  unit.  If 
Pdetect{x,y)  =  w,  we  know  that  if  a  target  stays  at  (x,y)  for 
l  time  units,  the  probability  of  being  detected  is  1  —  (1  —  wy. 
The  detection  probability  also  indicates  the  detection  delay. 
If  Pdetect  y)  =  w,  with  c%  confidence  level,  the  detection 
delay  is  within  ln(l  —  c%) /  ln(l  —  w)  time  units.  For  example, 
if  Pdetect{x,y )  =  0.5,  the  detection  delay  is  within  5  time 
units  with  95%  confidence  level.  If  the  system  requirement  is 
specified  by  the  desired  detection  performance,  the  designer 
can  compute  the  detection  probability  thresholds  at  different 
locations  based  on  the  relationship  between  the  detection 
probability  and  the  detection  performance. 

For  moving  targets,  we  are  interested  in  both  the  detection 
performance  and  the  tracking  performance.  Given  the  node 
detection  model  for  moving  targets,  the  relationship  between 
the  detection  performance  and  the  detection  probability  for 
moving  targets  is  the  same  as  that  for  static  targets.  The 
tracking  performance  mainly  depends  on  the  frequency  of  the 
detection  reports  from  the  network.  A  finer  grained  tracking  is 
achieved  if  the  time  interval  between  the  consecutive  detection 
reports  is  smaller.  The  frequency  of  the  detection  report  is 
indicated  by  the  detection  probability.  If  Pdetect{x ,  y)  =  w,  the 
average  time  interval  between  the  consecutive  detection  reports 
at  (x,  y)  is  \/w.  For  example,  if  the  detection  probability  in  an 
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area  is  0.5,  the  network  on  average  generates  detection  reports 
every  two  time  units  when  a  target  appears  in  that  area.  In  this 
case,  we  are  informed  the  locations  of  the  target  for  half  of  the 
time,  but  we  do  not  exactly  know  where  the  target  is  for  the 
other  half  of  the  time.  However,  if  the  detection  probability 
in  the  area  is  0.99,  we  get  the  detection  reports  almost  every 
time  unit.  In  this  case,  we  obtain  a  finer  grained  trajectory  of 
the  target,  together  with  the  duration  of  the  target  associated 
with  each  location  of  the  trajectory.  If  the  desired  tracking 
performance  at  different  locations  is  provided,  the  designer 
can  obtain  the  detection  probability  threshold  at  each  location 
based  on  the  relationship  between  the  detection  probability 
and  the  tracking  performance. 

V.  Problem  Formulation 

In  this  section,  we  formulate  the  node  deployment  problem 
as  an  integer  linear  programming  problem,  which  has  been 
proven  to  be  NP-hard.  Note  that  both  the  problem  formulation 
and  our  differentiated  deployment  algorithm  are  not  restricted 
to  the  detection  model  we  use.  While  we  use  the  detection 
model  described  by  Equation  (1),  any  detection  model  can 
be  applied.  We  first  assume  that  the  node  detection  model 
remains  the  same  no  matter  where  the  node  is  deployed.  A 
solution  is  provided  when  different  locations  result  in  different 
node  detection  models,  which  is  mainly  caused  by  terrain 
complexities,  such  as  obstacles  in  the  environment. 

In  this  paper,  we  study  the  system  that  transforms  the 
node  deployment  strategy  D  (the  input  of  the  system)  to 
the  logarithmic  collective  miss  probability  distribution  InM 
(the  output  of  the  system).  The  input-output  relationship  of 
the  system  can  be  characterized  by  the  node  detection  model 
p((m,n),(i,j)).  We  will  prove  later  in  this  section  that  this 
system  is  a  Linear  Shift  Invariant  (LSI)  system.  From  the 
property  of  the  LSI  system,  the  logarithmic  collective  miss 
probability  distribution  is  the  two-dimensional  convolution  of 
the  node  deployment  strategy  and  the  impulse  response  of  this 
system,  which  is  characterized  by  the  node  detection  model. 

In  order  to  simplify  the  relationship  among  the  node  deploy¬ 
ment  strategy  D,  the  node  detection  model  p((m,ri),(i,j)) 
and  the  collective  miss  probability  distribution  M,  we  first 
replace  p((m,n),(i,j))  with  pmisa((m,n),(i,j)).  As  shown 
in  Equation  (4),  the  collective  miss  probability  M(x,y )  is  the 
product  of  all  the  sensors’  probabilities  to  miss  the  target  at 
location  (x,y).  We  use  logarithm  to  convert  the  multiplicative 
relationship  shown  in  Equation  (4)  to  an  additive  relationship 
shown  in  Equation  (5). 

In  M(x,y)  =  Y  [D(i,j)xlnpmiss((x,y),(i,j))]  (5) 

(i,j)€_Grid 

For  simplicity,  we  use  I  to  denote  the  logarithmic  collective 
miss  probability  distribution  InM.  It  is  easy  to  see  that 
M (x,  y)  =  eI(x,y)  for  any  ( x ,  y)  £  Grid.  Through  the 
conversions,  the  system  described  in  Equation  (5)  is  indeed 
a  Linear  Shift  Invariant  System  (LSI),  which  takes  D  as  the 
input  and  I  as  the  output.  This  is  shown  by  the  following 
theorem. 


Theorem  1:  The  system  described  in  Equation  (5)  that 
transforms  the  node  deployment  strategy  D  to  I  is  a  Linear 
Shift  Invariant  (LSI)  System. 

Proof:  We  prove  that  the  system  described  in  Equation 
(5)  has  the  three  properties  of  an  LSI  system  as  follows: 

1.  Scaling.  Considering  two  inputs  1)\  and  D2,  and  D2  = 
XDi ,  in  which  A  is  an  arbitrary  nonnegative  integer  constant. 
Suppose  Ii  is  the  output  of  I)\  and  J2  is  the  output  of  D2. 
Then  for  any  ( x ,  y)  £  Grid ,  we  have: 

l2(x,y)=  Y  [-°2(i,i)  X  lnpmiss((x,!/),(i,i))] 

( i,j)€.Grid 

=A  Y  Pifhi)  x  j))] 

( i,j)GGrid 

=A/i  (x,  y) 

Therefore,  J2  =  \I\  and  the  system  conforms  to  the  scaling 
property. 

2.  Superposition.  Assume  that  the  outputs  of  two  arbitrary 
inputs  D\  and  D2  are  I\  and  /2,  correspondingly.  Suppose  /' 
is  the  output  D\  +  Z)2.  Then  for  any  (x,y)  £  Grid,  we  have: 

l'(x,y)=  Y  [(£*i(»>  j)  +  D2(i,  j))  X  lnpmiss((x,  y),  {i,  j))] 

Grid 

=  Y  [Dl  (iG)xlnpmiss((x,y),(i,j))] 

(i,j)€.Grid 

+  Y  XlnPmise(G,y),(i,j))] 

(i,j)€Grid 

y)  +  A(x,  y) 

Therefore,  /'  =  7)  +  /2  and  the  system  conforms  to  the 
superposition  property. 

3.  Shift  Invariant.  Assume  the  output  of  an  arbitrary  input  D 
is  I.  Shift  the  input  D  by  m  in  the  x-axis  and  n  in  the  y-axis 
to  form  a  new  input  D' ,  in  which  D'(x ,  y)  =  D(x  —  m,y  —  n ) 
and  m  and  n  are  arbitrary  integers.  Denote  the  corresponding 
output  of  D'  by  I’.  Then  for  any  (x,y)  £  Grid’,  in  which 
Grid'  is  the  sensor  field  after  we  shift  the  original  sensor 
field  Grid  by  m  in  the  x-axis  and  n  in  the  y-axis,  we  have: 

l\x,y)  =  Y  [D'G G')  x  lnpmi.At*.  y), 

(i'  ,j')GGrid' 

=  Y  [£>(»'  -  TO,  j'  -  n)  X  lnpmiaa((x,y),  (i\  j'))] 

(i'  ,j')&Grid' 

=  [D(i'  —  Tn,j'  —  n) 

(i'  ,j')GGrid' 

X  In  Pmiss  ((x  —  m,  y  —  n),  (i  —m,j  —  n))] 

=  Y  X  -  m,y  -  n),(i,j))\ 

( i,j)£Grid 

—  I (x  —  m,  y  —  n) 

Therefore,  the  system  conforms  to  the  shift  invariant  prop¬ 
erty. 

From  the  above  three  properties,  we  have  shown  that  the 
system  is  an  LSI  system.  ■ 

Since  it  is  an  LSI  system,  we  use  the  impulse  response  to 
characterize  how  an  input  D  is  transformed  into  the  output  I. 
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We  use  g  to  denote  the  impulse  response  in  this  LSI  sytem. 
When  we  use  the  node  detection  model  shown  in  Equation 
(1),  g(x,y)  is  expressed  in  the  following  equation: 

g(x,y)  =  (o,  o)) 

_  f  ln(l  -  e-a'/x2+«2)  if  y/ X2  +  y2  <  =  rs 
\  0  if  s/ x2  +  y2  >  re 

Further,  as  an  LSI  system,  the  output  I  can  be  represented 
by  the  convolution  result  of  the  input  D  and  the  impulse 
response  g ,  which  is  shown  in  Equation  (7). 

I{x,y)  =  D(x,y)  *  g(x,y) 

=  S  D(*>i)  x  v-j)  O) 

i  3 

Two  dimensional  convolution  is  a  well  studied  area  in 
image  processing.  We  apply  the  same  techniques  used  in 
image  processing  to  transform  the  convolution  to  the  ma¬ 
trix  multiplication  by  constructing  the  corresponding  matrices 
[18]  [19],  We  do  not  discuss  the  details  of  how  to  construct 
the  corresponding  matrices  in  this  paper  due  to  the  page  limit. 
Readers  are  referred  to  [  1 8] [19]  for  details.  By  using  the 
corresponding  matrices.  Equation  (7)  is  represented  by  the 
following  matrix  Equation: 

Ip  —  Gp  ■  Dp  (8) 

In  Equation  (8),  matrix  Gp  is  constructed  based  on  g,  matrix 
D.p  is  constructed  based  on  D  and  matrix  Ip  is  constructed 
based  on  I.  By  using  an  LSI  system  to  represent  the  rela¬ 
tionship  among  the  node  deployment  strategy  D,  the  node 
detection  model  p((m,  n ),  ( i ,  j ))  and  the  collective  miss  proba¬ 
bility  distribution  M,  we  finally  reach  a  simple  equation  which 
characterizes  the  relationship  among  all  three  elements.  Note 
that  if  TV  =  max(U ,  V),  both  the  dimensions  of  Ip  and  Dp  are 
TV2  x  1,  and  the  dimension  of  Gp  is  N2  x  TV2.  Also  note  that 
if  the  indices  in  Gp  start  from  zero,  the  elements  in  Gp  can  be 
interpreted  in  the  following  way:  the  value  of  Gp(x,y)  equals 
to  the  value  of  \npmiss({\x/N\,x%N),  ([y/N\,y%N)). 

By  now,  we  assume  that  all  the  nodes  share  the  same 
detection  model.  However,  this  assumption  does  not  hold  in 
certain  situations.  For  example,  the  sensing  capability  of  a 
node  may  be  undermined  by  the  obstacles  in  the  environment. 
This  problem  is  solved  by  modifying  the  matrix  Gp.  For 
example,  if  a  node  at  grid  point  (m,  n)  has  a  detection 
model  different  from  others,  we  change  the  values  of  the 
elements  in  the  (m  x  N  +  n)th  column  of  Gp.  If  the  value 
of  PmiSS({x,y),  ( m,n ))  is  changed  to  be  vnew,  Gp((x  x  N  + 
y),  ( m  x  N  +  n))  is  set  to  lnvnew.  In  this  way,  Equation  (8) 
allows  different  detection  models  at  different  locations  and 
the  impact  of  obstacles  is  incorporated.  By  allowing  different 
detection  models  at  different  locations,  the  system  is  no  longer 
an  LSI  system.  However,  it  is  still  a  linear  system  which 
can  be  proven  in  the  same  way  as  the  proof  of  the  scaling 
and  superposition  properties  of  Theorem  1.  So  Equation  (8) 


still  holds  and  both  our  problem  formulation  and  differentiated 
node  deployment  algorithm  are  not  affected. 

Based  on  Equation  (8),  we  are  able  to  formulate  the  node 
deployment  problem  as  an  integer  linear  programming  prob¬ 
lem.  We  use  uith  to  denote  the  miss  probability  threshold 
distribution  of  the  whole  field,  in  which  mth{x,y)  is  the  miss 
probability  threshold  at  location  (x,y).  Our  objective  is  to 
find  the  minimal  number  of  nodes  to  satisfy  that,  after  these 
nodes  are  deployed,  for  any  (x,  y)  £  Grid,  the  collective  miss 
probability  M(x,  y)  is  smaller  than  or  equal  to  mth{x,  y).  We 
set  I  to  In  mth-  Then  we  construct  the  corresponding  matrices 
for  D,  g  and  I.  The  problem  can  be  formulated  as  follows: 


Minimize:  sum(Dp) 

Subject  to:  Gp  ■  Dp  <=  I p 

The  elements  in  Dp  are  nonnegative  integers 


Fig.  1.  Integer  linear  programming  model. 

This  problem  is  a  typical  integer  linear  programming  prob¬ 
lem,  which  has  been  proven  to  be  NP-hard.  An  NP-hard 
problem  indicates  that  it  is  almost  impossible  to  find  the 
optimal  solution  when  the  input  size  is  large,  since  we  can 
not  afford  the  computation  time.  Therefore,  we  devise  an 
optimization  algorithm  called  DIET  DEPLOY  to  obtain  the 
suboptimal  result. 

VI.  Differentiated  Deployment  Algorithm 

In  this  section,  we  present  the  differentiated  deployment 
algorithm  DIFF_DEPLOY.  DIFF .DEPLOY  makes  use  of  the 
Linear  Shift  Invariant  (LSI)  property,  which  justifies  the  rep¬ 
resentation  of  the  relationship  among  the  node  deployment 
strategy,  the  node  detection  model  and  the  collective  miss 
probability  distribution  by  matrix  multiplication  as  shown  in 
Equation  (8).  Based  on  matrix  algebra,  it  is  clear  that  if  we 
know  the  miss  probability  threshold  distribution  and  the  node 
detection  model,  we  can  directly  compute  the  deployment 
strategy  Dp  based  on  Equation  (8)  as  follows: 

Dp  =  G~ 1  •  Ip  (9) 

The  result  Dp  computed  from  Equation  (9)  indicates  the 
number  of  sensors  we  should  place  at  each  location  to  satisfy 
that  M  =  mth »  based  on  the  node  detection  model.  If  we 
could  place  the  sensors  according  to  the  computed  result  Dp, 
the  total  number  of  nodes  deployed  would  be  minimal,  because 
Gp  ■  Dp  would  equal  to  Ip.  However,  the  computed  result  Dp 
does  not  reflect  the  real  world  situation.  In  reality,  the  number 
of  nodes  we  can  deploy  at  a  place  must  be  non-negative 
integers,  while  the  results  computed  from  Equation  (9)  can 
be  any  real  numbers,  including  negative  values.  Therefore,  we 
can  not  directly  use  the  result  Dp  as  the  deployment  strategy. 
However,  the  result  Dp  is  a  good  heuristics  to  decide  where 
to  deploy  sensors.  In  the  result  Dp,  a  larger  numeric  value 
requires  a  higher  number  of  nodes  to  be  deployed  to  satisfy 
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the  requirement.  Note  that  the  number  of  nodes  required  can  be 
any  real  number.  We  perceive  the  location  which  corresponds 
to  the  maximum  value  in  Dp  as  the  location  that  contributes 
the  most  in  satisfying  the  detection  requirement  if  a  sensor 
is  deployed  at  the  location.  When  we  decide  the  location  to 
deploy  the  next  sensor,  the  most  contributing  location  is  our 
first  choice. 


Procedure  DIFF_DEPLOY(Gp ,  Ip,  N){ 

UN  =  max(U,  V) 

invGp  =  inv(Gp)\  //Get  the  inverse  matrix  of  Gp 
remainlp  =  Ip] 

Dp  =  zeros(N  x  N,  1);  l*Dp  is  a  ( N  x  N)  x  1 

matrix,  initialized  with  all  zeros. 
while  sum(remainlp )  <  0  /*there  are  still  some 
locations  whose  miss  probability 
thresholds  are  not  satisfied*/ 
nextDp  =  invGp  x  remainlp ; 

Find  the  maximum  value  nextDp(k )  in 
nextDp,  which  satisfies  the  following 
constraints:  remainlp(k)  <  0 &&Dp(k)  ==  0; 
Dp{k)  =  Dp{k )  +  1;  /*Place  the  next  sensor 
at  location  k*l 

remainlp  =  Ip  —  Gp  x  Dp\  //Update  remainlp 
Set  any  positive  value  in  remainlp  to  zero; 
end 

} _ 


Fig.  2.  DIFF .DEPLOY  node  deployment  algorithm. 

DIFF.DEPLOY  makes  use  of  the  heuristics  provided  by  the 
result  Dp.  Figure  2  shows  the  pseudo  code  of  DIFF.DEPLOY. 
The  algorithm  runs  iteratively  and  places  one  sensor  in  the 
sensor  held  at  a  time.  During  each  iteration,  nextDp  is 
computed  based  on  the  remaining  requirement  remainlp, 
which  is  the  difference  between  Ip  specified  in  the  requirement 
and  the  logarithmic  collective  miss  probability  distribution 
provided  by  the  nodes  that  have  already  been  deployed.  The 
most  contributing  location  based  on  the  remaining  requirement 
remainlp  is  used  as  the  location  to  place  the  next  node, 
which  corresponds  to  the  maximum  value  in  nextDp.  Then 
the  remaining  requirement  remainlp  is  updated  based  on  the 
newly  deployed  node.  The  algorithm  terminates  when  all  the 
elements  in  remainlp  are  zero,  which  means  that  the  miss 
probability  thresholds  at  all  the  locations  are  satisfied. 

VII.  MIN_MISS  Algorithm  and  Uncertainty 

In  this  section,  we  briefly  describe  the  MIN  MISS  de¬ 
ployment  algorithm  presented  in  [10],  which  to  the  best 
of  our  knowledge  is  the  state-of-the-art  node  deployment 
algorithm  based  on  probabilistic  detection  model.  Although 
MINJVIISS  achieves  better  performance  than  the  random  node 
deployment  algorithm,  we  show  that  in  Section  VIII  our 
algorithm  DIM-  DEPLOY  achieves  much  better  performance 


than  MINJVIISS,  especially  for  differentiated  detection  re¬ 
quirements.  For  brevity,  we  use  RANDOM  to  denote  the 
random  node  deployment  algorithm. 

In  MINJVIISS,  for  each  iteration,  the  algorithm  decides 
the  location  for  the  next  node  to  be  deployed.  We  call  each 
location  at  which  no  sensor  node  is  deployed  a  “candidate 
location”.  If  we  deploy  a  new  node  at  candidate  location 
the  collective  miss  probability  at  location  (x,y)  is  then 
Pmiss{{x,y),{i,j))  x  M(x,  y).  The  overall  miss  probability 
Moverau(i,  j)  is  defined  as  the  sum  of  the  collective  miss 
probabilities  at  all  possible  locations  when  the  new  node  is 
deployed  at  location  or 

Moverau(i,j)  =  Pmiss{{x,y),(i,j))  X  M(x,y).  (10) 

( x,y)£Grid 

In  MINJVIISS,  candidate  location  (i,j),  which  has  the 
minimum  Moverau(i,  j)  among  all  the  candidate  locations, 
is  selected  for  the  next  sensor  to  deploy.  The  algorithm 
terminates  when  the  miss  probability  of  each  grid  point  is 
smaller  than  or  equal  to  its  miss  probability  threshold. 

The  paper  [10]  states  that  MINJVIISS  is  an  “uncertainty- 
aware  sensor  node  placement  algorithm”.  The  uncertainty  here 
means  that  when  we  deploy  sensors,  sensors  may  not  be  placed 
at  the  exact  intended  positions  because  of  uncontrollable 
conditions,  such  as  localization  error,  and  drifting  of  sensors  in 
undersea  sensor  networks  due  to  water  flows.  By  incorporating 
the  uncertainty  in  the  node  deployment,  the  node  detection 
model  is  modified.  However,  it  does  not  affect  the  deployment 
algorithm.  We  show  it  based  on  the  equations  from  [10]. 

In  [10],  a  Gaussian  probability  distribution  is  used  to 
model  the  deviation  between  the  intended  sensor  locations 
and  the  sensors’  actual  locations.  In  this  Gaussian  probability 
distribution  model,  the  intended  location  (x,  y)  is  the  mean 
value  and  the  standard  deviations  are  ax  in  the  x  dimension 
and  cjy  in  the  y  dimension.  With  this  model,  the  probability  of 
a  node  intended  to  be  deployed  at  location  (x,  y)  but  actually 
goes  to  location  ( x',  y ')  is: 

_  (x'-x)2  _  (y'-y)2 

2 a2  2^2 

.  .  g  x  y 

Pdepioy((x,y),(x  ,y  ))  =  - - -  (11) 

ZTUTx  Gy 

Let  Area  denote  all  the  possible  locations  a  sensor  can  be 
deployed  due  to  uncertainty.  Then  the  possibility  of  a  target  at 
grid  point  (m,  n)  being  detected  by  a  sensor  which  is  intended 
to  be  deployed  at  location  {x,  y)  can  be  denoted  in  Equation 
(12): 

Puncertainty{{m,n),  (x,y))  = 

E(x',,,')areaP((m'")-  G’ ,  V' ))Pdeploy  (G ,  V) ,  (x',v')) 
T,(x>y)£Area.Pdeploy((x,y),  (*',!/')) 

As  shown  in  Equation  (12),  by  introducing  uncertainty, 
the  node  detection  model  changes  from  p((m,n),(i,j))  to 
Puncertain({m,n),(i,j)).  We  call  the  new  model  that  con¬ 
siders  uncertainty  pUncertain{{m,n),(i,j)),  the  uncertainty 
aware  node  detection  model.  Figure  3  shows  the  uncertainty 
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aware  node  detection  model.  We  assume  that  ax  equals  ay, 
which  is  denoted  as  stdev  in  the  figure.  The  Gaussian  model 
is  truncated  and  we  set  the  maximum  distance  error  between 
the  intended  location  and  the  actual  location  to  3.  We  set  a 
to  0.5  to  compute  the  corresponding  detection  probabilities. 
As  shown  in  Figure  3,  a  larger  stdev  significantly  decreases 
the  detection  probability  near  the  intended  sensor  location,  but 
slightly  increases  the  detection  probability  far  away  from  the 
intended  sensor  location. 


Fig.  3.  Sensor  detection  model  when  uncertainty  is  considered. 

Since  MIN_MISS  is  a  general  node  deployment  algorithm 
as  DIFFJDEPLOY,  which  is  not  designed  just  for  one  spe¬ 
cific  node  detection  model,  it  makes  more  sense  to  decouple 
this  uncertainty  from  the  node  deployment  algorithm.  When 
uncertainty  is  introduced  into  the  node  deployment,  we  only 
need  to  feed  the  uncertainty  aware  node  detection  model, 
instead  of  the  original  node  detection  model,  into  the  node 
deployment  algorithm.  We  give  some  experimental  results  in 
Section  VIII-C  when  uncertainty  in  deployment  is  considered, 
to  show  the  robustness  of  our  proposed  algorithm. 

VIII.  Performance  Evaluation 

We  simulate  DIFFJDEPLOY,  MINJV1ISS,  and  RANDOM  in 
Matlab  and  compare  their  performance.  In  our  simulations,  we 
set  the  grid  dimension  to  50  x  50.  All  the  points  in  RANDOM 
are  the  mean  values  computed  from  20  trials.  Unless  otherwise 
specified,  the  sensing  range  is  set  to  7;  the  value  of  a  is  set 
to  0.5. 

For  this  performance  evaluation,  four  sets  of  experiments 
are  designed.  In  the  first  set,  we  compare  the  performance 
of  different  deployment  algorithms  for  uniform  detection  re¬ 
quirements.  In  the  second  set,  the  performance  of  different  de¬ 
ployment  algorithms  for  differentiated  detection  requirements 
is  evaluated.  In  the  third  set,  we  investigate  the  performance 
of  different  deployment  algorithms  when  uncertainty  in  the 
node  deployment  is  introduced.  In  the  forth  set,  we  compare 
the  performance  of  Dil  l-  DEPLOY  to  the  optimum  for  small 
scale  problems. 

A.  Evaluation  Set  1:  Uniform  Detection  Requirements 

We  first  study  the  performance  of  Dil  l-  DEPLOY, 
MINJVIISS  and  RANDOM  when  the  miss  probability  thresh¬ 
olds  at  all  locations  are  the  same. 

Figure  4  shows  the  performance  of  different  deployment 
algorithms  when  we  change  the  miss  probability  threshold. 


Fig.  4.  Performance  of  DIFFJDEPLOY,  MINJVIISS  and  RANDOM  with 
various  miss  probability  thresholds  for  uniform  detection  requirements. 

As  shown  in  Figure  4,  the  number  of  nodes  decreases  as  the 
miss  probability  threshold  increases  in  all  three  algorithms. 
Both  MINJVIISS  and  DIFFJDEPLOY  use  much  fewer  nodes 
than  RANDOM  for  uniform  detection  requirements.  However, 
DIFFJDEPLOY  achieves  better  performance  than  MIN_MISS 
even  for  uniform  detection  requirements.  In  DIFF .DEPLOY 
when  we  deploy  a  new  sensor,  both  the  remaining  miss 
probability  threshold  distribution  and  the  node  detection  model 
are  directly  integrated  into  the  computation  as  shown  in 
Equation  (9).  The  performance  advantage  of  DIFFJDEPLOY 
over  MINJVIISS  is  more  significant  when  the  miss  probability 
threshold  is  larger.  Overall,  DIFFJDEPLOY  uses  14%-42% 
fewer  sensors  than  MINJVIISS. 


Fig.  5.  Performance  of  DIFFJDEPLOY,  MINJVIISS  and  RANDOM  with 
various  a  values  for  uniform  detection  requirements. 

Figure  5  shows  the  performance  of  different  deployment 
algorithms  with  various  a  values,  when  the  miss  probability 
threshold  is  set  to  0.1.  As  shown  in  Figure  5,  the  number 
of  nodes  deployed  increases  as  the  value  of  a  increases, 
because  the  larger  the  value  of  a  is,  the  faster  the  detection 
probability  drops  when  the  distance  between  the  target  and 
the  node  increases.  Both  DIFFJDEPLOY  and  MIN  MISS 
greatly  outperforms  RANDOM  with  all  a  values.  Overall, 
DIFFJDEPLOY  achieves  20%-34%  better  performance  than 
MINJVIISS. 

Figure  6  shows  the  performance  of  different  deployment 
algorithms  when  we  change  the  sensing  range.  The  miss 
probability  threshold  is  set  to  0.1.  As  shown  in  the  figure, 
the  number  of  nodes  deployed  first  decreases  fast  in  all  three 
algorithms  as  the  sensing  range  increases.  Then  the  number 
of  nodes  starts  to  stabilize  when  the  sensing  range  is  above 
7,  because  the  detection  probability  becomes  small  when  the 
distance  between  the  target  and  the  sensor  is  large.  Overall, 


This  full  text  paper  was  peer  reviewed  at  the  direction  of  IEEE  Communications  Society  subject  matter  experts  for  publication  in  the  IEEE  SECON  2006  proceedings. 


Fig.  6.  Performance  of  DIFF_DEPLOY,  MIN_MISS  and  RANDOM  with 
various  sensing  ranges  for  uniform  detection  requirements. 


(a)  Without  Obstacle  (b)  With  Obstacle 

Fig.  7.  Miss  probability  threshold  distribution. 


DIFF_DEPLOY  uses  9%-34%  fewer  sensors  than  MIN_MISS 
with  various  sensing  ranges. 


the  issue  of  different  miss  probability  thresholds  at  different 
locations.  The  performance  advantage  of  DIFF_DEPLOY  over 
MINJVIISS  is  more  significant  for  differentiated  detection 
requirements.  By  adapting  to  the  different  miss  probability 
thresholds  at  different  locations,  Dil  l-  DEPLOY  uses  47%- 
57%  fewer  nodes  than  MINJVIISS. 


Fig.  9.  Performance  of  D IFF .DEPLOY,  MIN.MISS  and  RANDOM  with 
various  sensing  ranges  for  differentiated  detection  requirements. 

Figure  9  shows  the  performance  of  DIFF .DEPLOY, 
MINJVIISS  and  RANDOM  with  various  sensing  ranges. 
With  different  sensing  ranges,  the  performance  advantage  of 
DIFFJDEPLOY  over  MINJVIISS  is  also  much  more  significant 
for  differentiated  detection  requirements.  The  number  of  nodes 
used  in  DIFFJDEPLOY  is  50%-54%  fewer  than  that  used  in 
MINJVIISS. 


B.  Evaluation  Set  2:  Differentiated  Detection  Requirements 

We  study  the  performance  of  DIFF .DEPLOY,  MINJVIISS 
and  RANDOM  when  the  specified  miss  probability  thresholds 
at  different  locations  are  different  in  this  section.  Figure 
7(a)  shows  the  specification  of  the  miss  probability  threshold 
distribution  of  the  field.  Figure  8  and  Figure  9  show  the  results 
of  performance  evaluation  when  there  are  no  obstacles  in  the 
environment.  Then  we  show  their  performance  when  there  are 
obstacles  in  the  environment  in  Figure  10. 


Value  of  a 


Fig.  8.  Performance  of  DIFFJDEPLOY,  MIN.MISS  and  RANDOM  with 
various  a  values  for  differentiated  detection  requirements. 

Figure  8  shows  the  performance  of  the  three  deployment 
algorithms  with  various  a  values.  When  the  miss  probability 
thresholds  at  different  locations  are  different,  the  performance 
advantage  of  MINJVIISS  over  RANDOM  is  smaller  than  that 
when  the  miss  probability  thresholds  at  all  locations  are  the 
same,  because  MINJVIISS  does  not  take  into  consideration 


Value  of  a 


Fig.  10.  Performance  of  DIFFJDEPLOY,  MIN_MISS  and  RANDOM  with 
various  a  values  for  differentiated  detection  requirements. 

DIFF.DEPLOY  shows  promising  performance  for  differen¬ 
tiated  detection  requirements  when  there  are  no  obstacles.  We 
investigate  the  performance  of  DIFFJDEPLOY  when  there  are 
obstacles  in  the  environment.  The  distribution  of  the  obstacle 
is  shown  in  Figure  7(b).  We  assume  for  simplicity  that  a 
sensor  can  not  detect  a  target  if  there  is  an  obstacle  lying 
between  the  sensor  and  the  target.  Thus,  the  nodes  deployed 
near  the  obstacle  have  different  detection  models  than  others. 
We  show  the  performance  of  the  three  deployment  algorithms 
in  Figure  10  when  there  are  obstacles.  Compared  to  the 
performance  shown  in  Figure  8  when  there  are  no  obstacles, 
each  deployment  algorithm  uses  a  little  larger  number  of 
nodes.  The  existence  of  obstacle  restricts  the  sensing  ability  of 
the  nodes  deployed  close  to  it,  which  results  in  more  nodes. 
Overall,  DIFFJDEPLOY  requires  47%-54%  less  nodes  than 
MINJVIISS  with  various  a  values  when  there  are  obstacles. 
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C.  Evaluation  Set  3:  Uncertainty 

In  this  section,  we  show  the  performance  of  the  three  de¬ 
ployment  algorithms  when  uncertainty  in  the  node  deployment 
is  considered.  In  this  set  of  experiments,  we  use  the  same 
miss  probability  threshold  distribution  as  that  shown  in  Figure 
7.  Note  that  by  considering  uncertainty,  it  only  changes  the 
node  detection  model.  It  means  that  in  DIFFJDEPLOY,  only 
Yp  needs  to  be  replaced  with  the  one  that  is  constructed  based 
on  the  uncertainty  aware  node  detection  model. 
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Value  of  stdev 


Fig.  11.  Performance  of  DIFFJDEPLOY,  MIN.MISS  and  RANDOM  when 
uncertainty  is  considered. 


Figure  11  shows  the  performance  of  the  three  deployment 
algorithms  when  we  change  the  stdev  values.  In  this  set  of  ex¬ 
periments,  the  maximum  distance  error  between  the  intended 
location  and  the  actual  location  is  set  to  3.  With  uncertainty 
considered,  DIFFJDEPLOY  still  maintains  great  performance 
improvement  over  MINJVIISS.  Overall,  DIFF .DEPLOY  uses 
42%-54%  fewer  sensors  than  MIN  .MISS  with  various  values 
of  stdev,  when  uncertainty  is  considered. 


Maximum  Distance  Error 


Fig.  12.  Performance  of  DIFFJDEPLOY,  MIN_MISS  and  RANDOM  when 
uncertainty  is  considered. 


Figure  12  shows  the  performance  of  the  three  deployment 
algorithms  when  we  change  the  maximum  distance  error  be¬ 
tween  the  intended  location  and  the  actual  location.  The  stdev 
is  set  to  2  in  this  set  of  experiments.  As  shown  in  the  figure,  by 
varying  the  maximum  distance  error,  the  performance  of  each 
algorithm  does  not  change  much.  Overall,  DIFFJDEPLOY 
uses  45%-51%  fewer  sensors  than  MINJVIISS  with  different 
maximum  distance  errors. 

D.  Evaluation  Set  4:  Towards  Optimum 

We  have  shown  the  performance  superiority  of 
DIFFJDEPLOY  over  MINJVIISS  in  various  settings.  In 
this  section,  we  compare  the  performance  of  Dil  l-  DEPLOY 


to  the  optimum  for  small  scale  problems,  because  it  is 
computationally  infeasible  to  obtain  the  optimum  when  the 
input  size  is  large.  For  this  set  of  experiments,  the  miss 
probability  thresholds  at  different  locations  are  the  same;  the 
grid  dimension  is  set  to  5  x  5;  the  sensing  range  is  set  to  5; 
the  value  of  a  is  set  to  0.5.  We  use  OPTIMUM  to  represent 
the  optimal  node  deployment  algorithm,  which  always  uses 
the  minimum  number  of  nodes  to  satisfy  the  requirements. 


The  miss  Probability  threshold 


Fig.  13.  Performance  of  OPTIMUM,  D IFF .DEPLOY.  MINJVIISS  and 
RANDOM. 

Figure  13  shows  the  number  of  nodes  needed  for  OP¬ 
TIMUM,  DIFFJDEPLOY,  MIN  .MISS  and  RANDOM  when 
the  miss  probability  threshold  varies.  Overall,  for  small  scale 
problems,  DIFF .DEPLOY  is  the  algorithm  which  achieves  the 
closest  performance  to  the  OPTIMUM,  except  that  when  the 
miss  probability  threshold  is  0.01,  DIFFJDEPLOY  actually 
uses  one  more  sensor  than  MINJVIISS.  It  is  not  surprising 
to  see  that  in  some  special  cases,  DIFFJDEPLOY  uses  a 
few  more  sensors  than  MINJVIISS.  Both  DIFFJDEPLOY  and 
MINJVIISS  are  heuristic  based  algorithms  and  there  is  no 
guarantee  that  one  algorithm  always  performs  better  than  the 
other  one  under  all  settings.  For  this  small  scale  problem, 
DIFFJDEPLOY  shows  quite  close  performance  to  OPTIMUM, 
which  on  average  uses  1  more  sensor  than  OPTIMUM. 

IX.  Conclusions  and  Future  Work 

In  this  paper,  we  focus  on  differentiated  deployment  prob¬ 
lem,  in  which  the  required  detection  probability  thresholds  at 
different  locations  are  different.  We  apply  the  probabilistic 
detection  model,  which  is  more  realistic  than  the  binary 
detection  model.  We  show  that  the  relationship  between  the 
node  deployment  strategy  and  the  logarithmic  collective  miss 
probability  distribution  is  Linear  Shift  Invariant  (LSI).  We 
formulate  the  differentiated  deployment  problem  as  an  in¬ 
teger  linear  programming  problem,  which  is  a  well  known 
NP-hard  problem.  The  main  contribution  of  this  paper  is 
DIFFJDEPLOY,  a  differentiated  node  deployment  algorithm 
that  achieves  much  better  performance  than  the  state-of-the- 
art  algorithm  MINJVIISS,  especially  for  differentiated  detec¬ 
tion  requirements.  DIFFJDEPLOY  provides  the  flexibility  to 
enforce  the  differentiated  detection  requirements  and  at  the 
same  time  minimize  the  total  number  of  nodes  deployed. 

In  the  future,  we  will  investigate  other  applications  of  using 
the  LSI  relationship  between  the  node  deployment  strategy 
and  the  logarithmic  collective  miss  probability  distribution. 
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We  will  also  study  the  possible  solutions  when  the  sensing 
range  is  not  uniform,  in  which  case  the  sensing  ranges  in 
different  directions  are  different.  This  adds  extra  complexities 
into  the  current  problem,  which  requires  us  to  decide  not  only 
the  deployment  locations  of  the  sensors,  but  also  the  direction 
of  each  sensor. 
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